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I.  INTRODUCTION 

For  years  radar  developers  have  been  trying  to  evolve  computer  based  models  to  accurately 
represent  expected  radar  targets.  Traditionally,  targets  have  had  their  radar  “signature”  measured 
by  actually  recording  the  returned  signals  from  the  target  of  interest  and  building  a  database  for 
characterization  and  analysis.  This  is  an  expensive  process.  Each  individual  target  must  be  physi¬ 
cally  present  to  have  its  signature  characterized. 

Deriving  a  radar  signature  model  through  the  use  of  a  three-dimensional  target-drawing 
database  (not  requiring  the  physical  target  present)  has  long  been  the  dream  of  radar  developers. 
This  model  would  provide  a  cost  effective  alternative  to  real  data  collections  and  would  allow 
exploiting  the  full  capability  of  the  radar’s  operational  parameters.  Through  today’s  computa¬ 
tional  ability,  the  dream  of  an  automated  and  systematic  approach  to  theoretical  Radar  Cross 
Section  (RCS)  prediction  is  possible. 

This  paper  will  present  the  theories  behind  the  commonly  accepted  electromagnetic  tech¬ 
niques  used  for  RCS  prediction  and  enumerate  the  prevailing  computer  models  which  character¬ 
ize  RCS  for  high-frequency  complex  targets. 
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II.  RADAR  CROSS  SECTION  PREDICTION 


The  RCS  of  a  target  is  a  measure  of  power  reflected  back  toward  the  radar  when  illumi¬ 
nated  by  an  incident  wave.  RCS  gives  the  target  characteristics  without  including  the  effects  of 
transmitter  power,  receiver  sensitivity,  or  distance.  IEEE  defines  RCS  as  a  measure  of  reflective 
strength  of  a  target  defined  as  4k  times  the  ratio  of  the  power  per  unit  solid  angle  scattered  in  a 
specified  direction  to  the  power  per  unit  area  in  a  plane  wave  incident  on  the  scatterer  from  a 
specified  direction  [12].  The  formal  definition  is 

2 

o  =  lim  4.~iR2  ^sl 
*-»  N2 

where  R  is  the  distance  between  the  radar  and  target,  Es  is  the  field  strength  at  the  radar  receiver 
due  to  the  target  reflection  or  scattering,  and  E0  is  the  incident  field  strength  at  the  target  [11]. 
The  requirement  for  infinite  range  is  made  so  that  far  field  (plane  wave)  scattering  can  be  as¬ 
sured.  The  problem  of  predicting  RCS  becomes  one  of  determining  the  scattered  field  for  plane- 
wave  illumination. 

The  scattering  wave  breaks  into  three  regions:  (1)  the  low-frequency  or  Rayleigh  region 
where  the  wavelength  is  much  longer  than  the  scattering  body  size  and  the  scattering  process  is 
due  to  induced  dipole  moments  where  only  gross  size  and  shape  of  the  body  are  of  importance; 
(2)  the  resonant  region,  where  the  wavelength  is  on  the  same  order  as  the  body  size  and  the  scat¬ 
tering  process  is  due  to  surface  waves  and  optics;  and  (3)  the  high-frequency  region  where  the 
wavelength  is  much  smaller  than  the  body  and  the  scattering  process  is  principally  a  summation 
of  the  returns  from  isolated,  noninteracting  scattering  centers  [12].  These  regions  do  not  refer  to 
the  actual  frequency  of  the  incident  wave  but  to  the  size  of  the  target  when  compared  to  the  inci¬ 
dent  wavelength. 

The  calculation  of  the  RCS  for  a  complex  target  involves:  (1)  target  modeling  where  the 
geometry  of  the  target’s  exterior  surfaces,  ducts,  wedges,  apertures,  and  discontinuities  must  be 
described  in  terms  of  mathematical  models;  and  (2)  the  computation  of  scattered  electromagnetic 
fields  which  arise  from  different  mechanisms  such  as  specular  reflection,  diffraction,  and  multi¬ 
ple  scattering  [27].  To  model  the  above,  these  scattering  mechanisms  and  their  sources  must  be 
known. 
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HI.  RADAR  CROSS  SECTION  PREDICTION  TECHNIQUES 

The  basic  procedure  for  predicting  the  RCS  of  a  target  is  to  subdivide  the  target  into  a 
collection  of  simple  elements  such  as  flat  plates,  cylinders,  and  spheriods  and  assume  that  each 
scatterer  is  independent  of  all  the  other  scattered.  There  are  a  number  of  well  established 
theoretical  techniques  for  predicting  the  local  scattering  of  electromagnetic  radiation  for  high- 
frequency  RCS  prediction.  These  techniques  include  Geometrical  Optics  (GO),  the  Geometrical 
Theory  of  Diffraction  (GTD),  developments  from  GTD  such  as  the  Uniform  Theory  of  Diffrac¬ 
tion  (UTD),  Method  of  Equivalent  Currents  (MEC),  Physical  Optics  (PO),  and  the  Physical 
Theory  of  Diffraction  (PIT))  [15]. 

These  techniques  must  interact  with  one  another  to  be  able  to  predict  a  complete  RCS.  For 
example,  GO  fails  when  one  or  both  radii  of  curvature  become  infinite,  such  as  in  a  cylinder  or 
flat  plate,  so  we  would  have  to  use  PO  to  obtain  this  information.  The  techniques  cited  above  are 
addressed  in  this  section. 

A.  Geometrical  Optics 

The  GO  method  is  a  high-frequency  approximation  where  the  scattering  mechanism 
and  energy  propagation  are  accounted  for  by  classical  tubes  of  rays  [4].  It  is  best  used  for  ray 
tracing  that  does  not  encounter  shadow  boundaries  and  caustics.  Implementation  of  GO  requires 
that  there  exist  a  specular  point  on  the  body  and  that  the  principal  radii  of  curvature  be  stipu¬ 
lated.  The  principle  components  of  GO  include:  ray  paths,  ray  spreading,  and  reflection  coeffi¬ 
cients  [10]. 

GO  is  a  ray-tracing  procedure  where  the  wavelength  is  allowed  to  become  infinitesi¬ 
mally  small.  The  energy  propagates  along  slender  tubes  (rays)  according  to 

u  =  PeikS 

where  u  represents  the  complex  amplitude  of  a  magnetic  or  electric  field  component  transverse 
to  the  direction  of  propagation.  Magnitude  P  and  the  phase  factor  S  are  functions  of  position  in 
space  and  may  be  complex  numbers.  Propagation  is  in  the  direction  given  by  AS.  Surfaces  of 
constant  S  are  surfaces  of  constant  phase.  This  equation  is  not  valid  near  discontinuities  such  as 
edges  [11]. 

When  an  incident  ray  strikes  a  smooth  flat  surface  separating  two  media,  part  of  the 
energy  is  reflected  and  part  is  transmitted  across  to  the  second  medium.  Snell’s  law  indicates  that 
the  transmitted  ray  propagates  in  a  direction  different  from  that  of  the  incident  ray  (refraction), 
and  the  angle  of  the  reflected  ray,  as  measured  from  the  surface  normal,  is  equal  to  the  angle  of 
the  incident  ray  [11].  The  amplitude  and  phase  of  the  reflected  ray  are  given  by  Fresnel  reflec¬ 
tion  coefficients  [10]. 
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The  decay  or  increase  in  energy  can  be  calculated.  If  all  the  energy  which  entered  the 
tube  at  one  end  is  transmitted  to  the  other,  the  ratio  of  power  density  at  the  output  to  that  at  the 
input  is 


|Esl2  =  P1P2 
|E0|2  (s  +  Pi)  (s  +  p2) 

where  E(s)  is  the  field  intensity  at  the  output,  E(o)  is  the  field  intensity  at  the  input,  s  is  the  dis¬ 
tance  between  the  ends  of  the  ray  tube,  and  pi  and  p2  are  the  principal  radii  of  the  curvature  of 
the  output  wavefront  [11],  The  curvature  relationship  may  be  placed  in  the  ratio  of  power  densi¬ 
ty  at  the  output  to  that  at  the  input  and  then  into  the  RCS  formal  definition,  where  s  becomes  R. 
If  R  is  forced  to  infinity,  the  angular  dependence  on  the  local  angle  of  arrival  of  the  incident  ray 
disappears;  also,  the  angular  rotation  of  the  principal  planes  of  the  body  radii  of  curvature  also 
disappears.  The  result  is 

0  =  ;taia2 

where  aj  and  a2  are  the  principal  radii  of  the  body  at  the  specular  point.  One  of  the  problems  in 
computational  schemes  is  the  identification  of  the  specular  point. 

This  simple  equation  may  be  used  for  bistatic  as  well  as  monostatic  cases  as  long  as 
there  is  no  forward  direction.  It  may  be  used  for  metallic  bodies,  but  the  specular  reflection  must 
be  multiplied  by  the  reflection  coefficient  if  the  body  is  dielectric.  Dielectric  bodies  are  charac¬ 
terized  by  rays  that  penetrate  the  body  as  well  as  those  reflected  by  the  surface  [11]. 

In  GO  approximation,  the  region  of  space  is  divided  into  distinct  illuminated  and  shad¬ 
ow  regions;  therefore,  this  method  cannot  account  for  any  diffraction  effect  [4].  Also,  GO  cannot 
take  care  of  the  phase  of  the  field,  caustic  regions,  shadow  boundaries,  or  shadow  region  fields 

[4]- 


B.  Geometrical  Theory  of  Diffraction 

J.  B.  Keller  extended  GO  in  GTD  to  include  regions  where  diffracted  fields  are  impor¬ 
tant,  such  as  in  shadowing  [4].  GTD  contains  diffracted  ray  paths,  including  surface  rays,  and 
diffraction  coefficients  [10].  Figure  1  displays  examples  of  diffracted  rays  for  edge,  wedge, 
comer  tip,  multiple  diffracted,  and  curvature  [10], 
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/  Round  bodies  Tip  or  apex 

Figure  1.  Diffraction  Mechanisms 

Reference  10  shows  that  the  general  form  used  to  implement  GTD  is  to  express  the 
scattered  field  as  a  product  of 

E  =  T  Ad(s)eiksn-  E 

where 

T  =  Reflection  coefficient, 

Ad(s)  =  Ray  spreading  factor, 

eiks  _  phase  variation  relative  to  the  specular  point, 

D  =  Diffraction  matrix  which  relates  the  incident  to  the 
scattered  electric  field  at  the  diffraction  point. 


where 

Du,  Dj.  =  Scalar  diffraction  coefficients. 
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Shown  in  Reference  4,  the  D  is  based  on  exact  solution  of  the  canonical  problem  of  diffraction 
such  as  the  wedge  diffraction  coefficients: 


where 


Ds,h  ~ 


(2/n)  sin(:r/n) 
#M)  sin  Po 


cos 


n  ) 


_ 

_ 

-r 

(K 

\  (<b  +  <M 

cos\n 

)  005  In) 

> 

- 

\  / 

n  =  (2k  -  a)/K  with  a  as  the  interior  wedge  angle. 

p0  =  The  smaller  angle  between  the  direction  of  incidence 

and  the  tangent  at  the  point  of  incidence. 

<j)0,  <|>  =  The  angles  of  incidence  and  scattering,  respectively, 

s,  h  =  Soft  and  hard  wedges;  s  and  h  correspond  to  plus 

and  minus,  respectively. 


Keller’s  diffraction  coefficients  [4]  for  a  straight  wedge  are  shown  in  Figure  2. 


Figure  2.  Diffraction  of  a  Plane  Wave  by  an  Infinite  Plane  Wedge 

Several  limitations  have  been  identified  using  GTD.  These  include:  (1)  the  diffracted 
field  becomes  infinite  at  both  incident  and  reflected  shadow  boundaries;  (2)  the  diffraction  coef¬ 
ficient  becomes  infinite  at  the  edge,  thereby  breaking  the  edge  condition;  (3)  the  expression  for 
diffraction  coefficient  has  been  acquired  only  by  comparison  with  an  exact  solution  of  a  canoni¬ 
cal  problem  and  does  not  follow  as  an  integral  part;  (4)  the  higher-order  terms  in  the  expression 
for  a  diffracted  wave  cannot  be  determined;  and  (5)  the  solution  becomes  singular  at  a  caustic 

[4]- 


6 


C.  Uniform  Theories  of  Diffraction 

Two  different  theories  have  been  developed  to  handle  the  shadowing  limitation  ob¬ 
served  in  GTD.  Reference  1  shows  the  development  of  the  Uniform  Asymptotic  Theory  (UAT) 
of  diffraction.  This  was  on  the  assumption  that  the  field  solution  in  an  edge  diffraction  problem 
can  be  expanded  in  a  particular  asymptotic  series  implicating  a  Fresnel  integral.  The  UAT  solu¬ 
tion  is  uniformly  valid  close  to  and  away  from  the  shadow  boundaries  [4]. 

The  most  popular  method  is  UTD  which  was  developed  by  Kouyoumjian  and  Pathak 
[4].  UTD  can  handle  large-sized  scattering  objects  in  an  efficient  manner.  The  larger  the  object 
in  terms  of  a  wavelength,  the  better  the  solution.  Complex  objects  can  be  systematically  built  up 
using  the  major  local  scattering  centers  based  on  diffraction  coefficients  [22].  The  singularities 
in  the  diffraction  coefficients  are  defeated  by  multiplying  the  diffraction  coefficients  by  a  Fresnel 
integral.  At  the  shadow  boundary,  or  the  reflection  boundary,  the  modifying  Fresnel  integral  is 
zero  while  the  diffraction  coefficient  is  infinite.  The  product  of  the  two  stays  finite  [10]. 

Figure  3  shows  the  diffraction  of  an  electromagnetic  wave  by  a  curved  wedge  [4]. 


A 

n 


Figure  3.  Diffraction  by  a  Curved  Wedge 

Reference  4  gives  the  diffraction  coefficients  uniform  at  and  around  the  shadow  boundaries  as: 
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Ds’h(4>,  4)',  p0) 


—  exp(i:t/4) 


2n  jlxk  sin  po 
F[kLa+  (4>  -  4>')]  cot  — 

+  F[kLa-  (4,  -  4,')]  cot  — 


F[kLa+  (4,  +  4>')j  cot 
+  F[kLa“  (4>  +  4>')]  cot  — — ^ — — 


where 


F(x)  =  2i  Jx  exp1 


f 


exp  1x2  dt 


2  2nJtN±  -  (p) 


a±(P)  =  2  cos 


where  N*  are  the  integers  that  most  nearly  satisfy  these  equations: 

2rcnN+  -  (P)  =  Jt  with  p  =  +  <(>', 

2jtnN--(P)  =  k  with  p  =  4>  —  4'. 

L  is  a  distance  parameter  defined  as 

L  =  S2  sin  p0  for  plane  wave  incidence, 


L 

L 

r' 

r 

S 

N± 

Po 


^  ,  for  cylindrical  wave  incidence, 


ss 


s  +  s 


7  sin2  po  for  conical  and  spherical  wave  incidence, 


radius  of  curvature  of  the  cylindrical  wave  normally  incident  on  the  edge, 
perpendicular  distance  of  the  field  point  from  the  edge, 
distance  along  the  ray  from  the  edge  to  the  observation  point, 
functions  of  p  and  n, 

the  smaller  angle  between  the  direction  of  the  incidence, 


and  the  tangent  is  drawn  at  the  point  of  incidence  on  the  edge  [4]. 

F(x)  is  a  modified  Fresnel  integral  which  is  shown  approximated  as  [4]: 

F(x)  =  Jkx  -  2x  exp(irt/4)  -  j  x2  exp(-  ijt/4)J  exp^2  +  xj 


F(x)  =  ( 


1  + 


J_  _  lx  _  ;15  J_  ,  75 
2x  4x2  8  x2  16  x4 


for  x  <  1 


for  x  >  1 . 
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When  the  source  is  near  the  edge,  the  diffraction  coefficient  needs  a  correction  consid¬ 
ering  the  variation  of  incident  field  intensity  along  the  edge  of  the  wedge.  This  is  known  as  the 
uniform  slope  diffraction  correction.  Reference  4  defines  the  uniform  slope  diffraction  coeffi¬ 
cient  as: 

dDs.hlfo  4>';  go)  =  _  e~*/4 

d4>'  4n2  Jink  sin  (30 

x  |csc2  Jt  +  —  2D^'~  Fs[kLa+(4)  —  4>,)j 
_  csC2  "  ~  ~  *')  Fs[kLa-  (4.  -  4>’)]  - 

J  - 

±  ■  CSC2  *  +  Fs[kLa+  (<(>  +  <>')] 

-  CSC2  — — ^D+  —  1  Fs[kLa"  (<?  +  <j>')]  ■ 

J  J  + 

where 

Fs(x)  =  2ix  +  4x3/2eix  e_k2  dt 
or 

Fs(x)  =  2ix[l  -  F(x)]. 

At  a  shadow  boundary,  the  following  limit  is  used: 

CSC2  Fs[kl-a±®]  =  4ikLn2  |l  -  £  e-/4j 

where 

P  =  2jrnN±T(Jt  —  e) 
and  £  is  small. 

The  disadvantages  of  UTD  hinge  on  the  fact  that  it  is  assumed  that  the  most  important 
mechanisms  are  known.  If  some  terms  are  left  out,  the  accuracy  of  the  total  solution  may  di¬ 
verge;  accuracy  also  depends  on  knowing  the  diffraction  coefficients.  As  more  terms  are  added 
to  the  solution,  the  more  time  it  takes  to  trace  and  shadow  the  rays  [22].  UTD  cannot  handle 
caustics  nor  the  treatment  of  cross  polarization  [4]. 
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D.  Method  of  Equivalent  Currents 

MEC  is  an  analytical  tool  using  fictitious  currents  to  overcome  the  caustic  difficulties. 
MEC  has  the  ability  to  predict  the  scattered  fields  in  the  direction  of  a  caustic,  to  calculate  the 
diffracted  fields  from  an  edge  of  finite  length  and  for  observation  angles  away  from  the  cone  of 
diffracted  rays,  and  to  find  the  reflected  fields  from  curved  bodies  of  finite  length  [19]. 

The  basic  approach  is  to  assume  the  existence  of  electric  and  magnetic  currents  at  each 
point  around  the  edge  contour  and  sum  them  in  the  far  field  radiation  integral  [12] 


Ed 

% 

A 


-  ik'P, 


■ 

Z0leS  X  ^S  X  tj  +  lm^S  X  tjj 


,  —  ikr-s 


dt 


Far  field  Green’s  function, 

Unit  vector  aligned  along  the  contour. 


Millar,  Ryan  and  Peters,  Knott  and  Senior,  and  Michaeli  have  used  equivalent  currents 
in  many  of  their  studies.  Equivalent  currents  proposed  by  Ryan  and  Peters  are 


le  =  i2(X  -  Y)  (t-Ejj/kZo 
lm  =  i2(X  +  Y)  (t-Hjj/kZo  . 

where 


(1/n)  sin  (x/ n) 
cos  (Jt/n)  -  cos[(\Fs  -  Wj)/n 


Y  =  (1/n)  sin  (jt/n) 

cos  ( jt/n )  -  cos^s  +  Wj)/n] 

n  =  Exterior  wedge  angle  normalized  with  respect  to  x. 

'P;,  Ws  =  Angles  of  the  transverse  components  of  the  incident 
and  diffracted  directions  with  respect  to  one  of  the 
surfaces  meeting  at  the  edge. 

When  these  currents  are  used  in  the  far  field  radiation  integral,  all  scattered  fields  are  obtained. 
To  extend  MEC  to  general  bistatic  cases,  Knott  and  Senior  interpreted  the  product: 

sin2  p  =  sin  Pi  sin  ps 

where  i  and  s  denote  the  angle  subtended  by  the  incident  and  scattering  directions.  Reference  12 
denotes  when  the  equivalent  currents  are  placed  in  the  far  field  radiation  integral,  the  result  is  a 
contour  integral: 


Ed  = 


-  2EC 


f  eikU;-;) 

°^°Jcsin  pi  sin  ps 
(t-ej)  (X  -  Y)  s  x  (s  x  t)  +  (t-hj)  (X  +  Y)s  x  tj  dt 


e;  =  Unit  vector  aligned  along  incident  electric  field. 

A 

hj  =  Unit  vector  aligned  along  incident  magnetic  field. 
C  =  Illuminated  portions  of  the  edge. 
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Michaeli  devised  a  completely  different  set  of  diffraction  coefficients  considering  the 
far  field  contribution  due  to  two  narrow  surface  strips  (one  on  each  face  of  a  wedge)  meeting  at 
an  edge  element  and  integrated  the  induced  currents  along  the  strips. 

A  disadvantage  of  MEC  is  the  complexity  of  the  diffraction  coefficients.  This  may 
prohibit  the  implementation  of  the  method  for  complex  bodies  which  is  composed  of  hundreds 
of  facets  [4]. 


E.  Physical  Optics 

PO  overcomes  the  infinities  of  flat  and  singly  curved  surfaces  prominent  in  GO  by 
approximating  the  induced  surface  fields  and  integrating  them  to  obtain  the  scattered  field  [12]. 
PO  is  a  high-frequency  technique  where  the  target  is  replaced  by  induced  surface  currents  that 
serve  as  the  source  of  scattered  fields.  In  obtaining  the  currents,  it  is  assumed  that:  (1)  the  radii 
of  curvature  of  the  surface  are  large  compared  to  wavelength,  and  (2)  currents  exist  only  in  the 
area  that  is  directly  illuminated  by  the  incident  wave  and  currents  on  the  illuminated  surface 
have  the  same  characteristics  as  those  on  the  infinite  plane  tangent  to  the  surface  at  the  point  of 
incidence  [4].  PO  relies  on  the  estimates  of  the  surface  fields  from  GO  [10]. 

Reference  14  noted  twelve  different  interpretations  of  PO.  These  depend  on  the  choice 
of  current  distribution  on  the  target  surface  and  other  characteristics.  Reference  4  gives  three 
possible  choices  of  the  currents  on  the  target  surface  as: 

E  -  H-Formulation:  Js  =  n  X  ^H1  +  Hs);  M  =  (e1  +  Es)  x  n 

H-Formulation:  Js  =  2^n  X  H‘);  M  =  0 

E-Formulation:  Js  =  0;  M  =  2^ES  x  n) 


where  E\  H’  and  Es,  Hs  are  the  incident  and  scattered  fields  at  the  target  surface  and  are  related 
by  boundary  conditions. 

In  general,  PO  current  is  imbedded  as  an  unknown  in  an  integral  equation  whose  solu¬ 
tion  provides  the  required  current  for  field  calculations.  Reference  4  explained  that  ignoring  the 
effects  of  any  surface  current  in  the  direction  of  observation,  the  Stratton-Chu  equations  are 
modified  to  give  the  scattered  fields: 


Es 


Hs 


+  jL_koI  I  s  x  (m  +  T)0s  x  Js)  expj^ikor'-  (i  -  sj 

J  S 

_  f  £  x  (js  +  ^s  x  m)  expj^  -  ikor'-  (i  - 


ds' 


where 


i,  s  =  Unit  vectors  in  the  directions  of  incidence  and  observation. 
S  =  Illuminated  portion  of  target  surface. 

All  other  parameters  are  shown  in  Figure  4  [4]. 


ll 


2 


P(r,  0,  0) 


Figure  4.  Electromagnetic  Scattering  from  a  Surface 

The  two  equations  above  may  be  considered  the  generalized  PO  formulation  and  are  used  to 
obtain  the  scattered  field  with  a  knowledge  of  the  currents  Js  and  M  obtained  either  by  solving 
appropriate  integral  equations  or  by  suitable  approximations  [4]. 

The  accuracy  of  PO  can  be  improved  by  improving  the  assumed  current  distributions. 
In  PO,  the  contributions  from  edges  are  ignored.  Disadvantages  include:  (1)  cannot  predict  cur¬ 
rents  due  to  a  discontinuity  on  a  body;  (2)  cannot  handle  depolarization  for  a  monostatic  single 
scatterer;  (3)  cannot  predict  diffraction  effects  (shadow  region  fields);  and  (4)  is  considered  inac¬ 
curate  for  wide-angle  lobe  predictions  [4]. 

F.  Physical  Theory  of  Diffraction 

A  more  accurate  representation  of  the  scattered  fields  than  PO  was  developed  by 
Ufimtsev,  known  as  PTD.  PTD  was  developed  for  treating  edges.  Ufimtsev  represented  the 
scattered  field  as  the  sum  of  the  PO  contribution  and  an  edge  contribution,  using  the  exact 
solution  of  the  two-dimensional  wedge  problem  [12]. 

Using  Sommerfeld’s  solution  of  the  wedge  problem,  Ufimtsev  represented  the 
incident  field  plus  the  scattered  fields  as: 

Ez  =  Eoz  [u(r,  VTS  -  W,)  -  u(r,  Ws  +  Wj] 

Hz  =  Hoz  [u(r,  Ws  -  W,)  +  u(r,  Ws  +  ¥•)] 

where 


u(r,  W] 


JL  (  _ g  —  ikr  cos  p  _  d 

2a  Jc  1  -  exp  irt((3  +  W)/a] 


where  r  is  the  distance  from  the  edge  to  the  point  of  observation,  'Ps  is  the  angular  coordinate  of 
that  point  above  one  face  of  the  wedge,  T7;  is  the  direction  of  arrival  of  the  incident  wave,  a  is 
the  external  wedge  angle,  and  C  is  the  Sommerfeld  contour  in  the  complex  plane  [12]. 
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Depending  on  the  direction  of  arrival  of  the  incident  wave,  there  are  three  possible  PO 
contributions  to  subtract  from  the  exact  solution.  One  occurs  if  the  upper  face  is  illuminated,  but 
not  the  lower  face;  another  occurs  if  the  lower  face  is  illuminated,  but  not  the  upper;  and  the 
third  occurs  if  both  faces  are  illuminated  [12].  Results  will  be  given  in  different  forms  for  the 
three  different  contributions.  Ufimtsev’s  results  are  summarized  below: 


=  Ef 


ei(kr+j[/4) 

/2jrkr 


tt  „  e^  +  rc/4) 

02  y2jtkr 


where  the  diffraction  coefficients  are 


(X  -  Y)  -  (Xx  —  Yj)  0  <  Wj  <  a  —  n 

(X  -  Y)  -  (Xj  -  Yj)  -  (X2  -  Y2)  a  -  n  <  ^  <  n 

(X  -  Y)  -  (X2  -Y£  n  <  ^  <  a 


'(X  +  Y)  -  (Xt  +  Yj)  0  <  <  a  -  n 

g  =  .  (X  +  Y)  -  (Xx  +  Yj)  -  (X2  +  Y2)  a  —  n  <  ^  <  it 

(X  +  Y)  -  (X2  +  Yj)  n  <  Wi  <  a . 

Reference  12  describes  the  unsubscripted  coefficients  as  Keller  coefficients  shown  on  page  13  of 
this  report.  The  subscripted  coefficients  are  due  to  the  uniform  PO  contributions. 

X,  =  -  1  tan  [(V,  -  V,)/2] 

Y,  =  -  i  tan  [(V,  +  V,)/2] 

X2  =  \  tan  [(¥.  -  VO/2] 

V2  =  -  \  tan  [a  -  {V,  +  Vj)/2] 

A  disadvantage  in  using  PTD  is  that  it  yields  only  the  contributions  due  to  edges. 

Some  other  method,  such  as  PO,  must  be  used  to  obtain  the  surface  contributions  [12]. 
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IV.  CURRENT  MODELING  APPLICATIONS 


Distinct  computer  codes  exist  today  for  high-frequency  RCS  prediction.  A  major  issue  for 
code  developers  has  been  in  the  implementation  of  electromagnetic  theory  for  practical  applica¬ 
tions  [13].  The  target  is  essentially  divided  into  two  parts:  how  the  incident  and  scattered  fields 
are  related  at  a  particular  point  on  the  target,  and  the  interaction  and  combination  of  the  scat¬ 
tered-field  contributions  from  different  points  of  the  target  [15]. 

Code  development  began  in  the  early  1970’s.  One  such  code  was  MISCAT,  developed  by 
Northrop  Corporation  [18].  High-frequency  computer  codes  have  been  ameliorated  over  the  past 
several  years.  The  most  prevailing  codes  today  include:  XPATCH,  RCS  -  Basic  Scattering  Code 
(RCS-BSC),  Simulated  Radar  Image  (SRIM),  and  MAX-TRACK. 

A.  XPATCH 

XPATCH  is  a  general  purpose  RCS  prediction  code  used  to  calculate  the  polarimetric 
radar  return  from  different  targets.  The  code  was  developed  by  DEMACO,  Inc.  under  the  joint 
tri-Service  sponsorship  of  Wright  Laboratory  at  Wright-Patterson  AFB,  Ohio,  Phillips  Labora¬ 
tory  at  Kirtland  AFB,  New  Mexico,  the  Army  Research  Laboratory  at  Ft.  Belvoir,  Virginia,  and 
the  Naval  Air  Warfare  Center  at  Pt.  Mugu,  California  [2].  The  first  version  was  written  in  1988; 
however,  it  wasn’t  until  late  1991,  after  further  development,  that  it  was  used  in  RCS  prediction 
[13]. 

XPATCH  consists  of  three  parts:  electromagnetics  which  includes  XPATCH1, 
XPATCH2,  XPATCH3,  and  XPATCH4;  CAD  and  visualization  tools;  and  a  graphical  user  inter¬ 
face  [13].  All  codes  within  XPATCH  may  be  used  as  independent  units.  ACAD  is  the  preferred 
CAD  system  to  use  for  developing  input  target  files;  however,  other  CAD  systems  may  be  used 
with  XPATCH.  There  is  a  translator  entitled  “FRED”  which  will  convert  from  solid  geometry  to 
surface  geometry  (facets)  or  the  reverse. 

XPATCH1  uses  triangular  facet  files  for  input  target  geometry.  While  this  format  is 
common  for  air  targets,  many  existing  CAD  models  for  ground  vehicles  were  built  using  BRL- 
CAD.  BRL-CAD  uses  combinatorial  solid  geometry  where  different  geometric  primitives 
(boxes,  cones,  ellipsoids,  etc.)  are  combined  using  Boolean  operations  such  as  union,  intersec¬ 
tion,  and  subtraction.  XPATCH2  uses  the  solid  geometry  models  [13]. 

SAR  images  were  being  developed  using  the  traditional  FFT  method.  For  each  SAR 
image,  the  scattered  field  from  the  target  was  typically  calculated  for  256  incident  angles  and 
256  frequencies  before  the  2D  Fast  Fourier  Transformation  (FFT)  was  taken.  Computer  time  was 
extensive.  A  direct  shooting  and  bouncing  ray  technique  was  implemented  in  XPATCH3,  which 
uses  triangular  facet  files  for  target  geometry,  and  XPATCH4,  which  uses  the  solid  geometry 
target  files  as  input  [13]. 

XPATCH  is  currently  based  on  the  Shooting  and  Bouncing  Ray  (SBR)  technique. 
Reference  3  has  recognized  the  advantage  of  SBR  and  have  pursued  establishing  its  validity.  In 
SBR,  a  dense  grid  of  rays  is  emitted  from  the  radar  direction  toward  the  target.  These  rays  are 
traced  according  to  the  GO  theory  as  they  bounce  from  one  part  of  the  geometry  to  another  and 
then  exist.  The  tracing  includes  the  effects  of  polarization,  ray  divergence  factor,  and  layered 
material  reflection.  At  the  point  where  a  ray  exits  the  target,  a  frequency  domain  or  time  domain 
PO  integration  is  done  to  calculate  the  scattered  far  field  from  the  target.  All  single  and  multiple 
bounce  contributions  are  accounted  for  by  the  GO  or  PO  theory  [13].  This  ray  bouncing  tech¬ 
nique  is  completely  adequate  for  an  aperture  opening  large  compared  to  the  wavelength;  howev¬ 
er,  as  the  opening  is  decreased,  additional  features  of  the  actual  field  must  be  accounted  for  in 
Reference  16. 
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Reference  16  indicates  the  bouncing  technique  that  is  carried  out  in  three  parts: 

1.  Given  the  geometry  of  the  cavity  and  the  incident  field,  the  ray  paths  in  the 
cavity  must  be  found  by  ray  tracking. 

2.  Determine  the  field  amplitude  of  the  exit  rays  on  the  aperture  based  on 
geometrical  optics.  The  ray  tube  divergence  factors  and  the  reflection 
coefficients  must  be  calculated. 

3.  Determine  the  backscattered  field  and  the  RCS  using  Kirchhoff’s  approximation 
(physical  optics). 

Following  is  a  general  description  of  these  three  steps. 

Parallel  rays  are  launched  from  the  incident  direction.  Each  ray  is  represented  by  a 
line  in  space  with  a  reference  point  (x0,  y0,  Zo)  and  a  direction  vector  (sx,  Sy,  sz).  Any  point 
(xi,  yi,  zi)  along  this  line  is  described  by 

(xi,  yi,  zi)  =  (x0,  y0>  Zo)  +  (sx,  Sy,  sz)t 

where  t  is  time.  If  the  phase  of  the  field  at  (xj,  yi,  z^)  lags  point  (xo,  y0,  z0),  t  is  a  positive 
quantity.  The  direction  vector  of  the  incident  rays  is  given  by 

sx  =  -sin  0*  cos  <j)' 

sy  =  -sin  01  cos  4*1 

sz  =  -cos  O'. 

The  reference  point  (x0,  y0,  z0)  on  the  incident  plane  21  can  be  related  to  the  point  (xa,  ya,  0)  on 
the  aperture  2a  by 

x0  =  (sy  +  S 2)xa  -  sxsyya  +  sxtD 

y0  =  (si  +  s£)ya  -  sxsyxa  +  syt0 

v 

z0  =  Sxszxa  SySzya  +  szt0 . 

The  t0  above  determines  how  far  21  is  from  the  aperture  plane  and  can  be  chosen  arbitrarily  [10]. 

Once  the  incident  rays  have  been  defined,  the  impact  point  of  each  ray  on  the  inner 
wall  of  the  cavity  can  be  determined.  This  is  achieved  by  the  equation  describing  the  cavity 
z  =  (f(s,  y))  and  any  point.  For  example,  if  the  cavity  is  a  circular  cylinder  with  radius  a,  the 
intersection  is  found  by  substituting  into  the  equation  for  any  point  along  the  line 

t  =  (-B  +  7b2  -  4Acj/2A 

where 

A  =  sx  +  s^ 

B  =  2(sxx0  +  syy0) 

C  =  Xq  +  yo  ~  a2  . 
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The  reflected  ray  must  satisfy  Snell’s  law,  such  as  it  must  lie  in  the  plane  of  incidence 
and  the  angle  of  reflection  must  equal  the  angle  of  incidence.  From  Figure  5,  Reference  16  de¬ 
fines  a  unit  vector  by 


m  =  |l0  x  nj/ sin  0c 
m  =  Perpendicular  to  the  plane  of  incidence. 

A 

10  =  Unit  vector  pointing  from  point  1  to  point  .0. 

The  departure  of  the  ray  from  the  cavity  is  easily  detected  when  the  ray  intersects  the  aperture 
2a-  For  every  ray  launched,  a  set  of  impact  points  inside  the  cavity  and  the  direction  of  the 
exiting  ray  are  obtained  [16]. 

Once  the  ray  paths  inside  the  cavity  are  found,  the  field  amplitude  can  be  established. 
For  each  ray  launched,  a  set  of  impact  points  is  found  inside  the  cavity'.  The  aperture  field 
associated  with  the  exit  ray  can  be  determined.  Using  GO,  the  electric  field  obeys  the  following 
relationship: 


where 


E  (xI+i»  yj+r  zI+i) 


(DF)i-(r).-E(xI,  y?  z-)*e  —  j<Phase) 


phase 


+  (yi+i 


E  fxj,  yT,  Zjj  =  Incident  Field  at  (x*,  yj,  Zj). 


(DF)i 


Planar  reflection  coefficient  matrix  at  the  ith  reflection 
point  where  the  original  curved  interface  (Fig.  5)  is 
replaced  by  its  tangent  plane  at  the  reflection  point. 

Divergence  factor  which  governs  the  spreading  of  the 
differential  ray  tube  from  just  after  the  ith  reflection 
to  just  before  the  (i  +  1)*  reflection. 
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A 


Figure  5.  Local  Coordinate  System  of  a  First  Impact  Point 

The  remaining  task  for  amplitude  tracking  is  to  determine  the  planar  reflection  coefficients  and 
the  incident  field. 

In  determining  the  backscattered  field  and  the  RCS,  PO  is  used.  The  outgoing  field  is 
replaced  by  an  equivalent  magnetic  current  sheet 

2E  (xN,  yN,  0)  x  z;  over  2A 
0  outside  2A 


Ks  =  Radiates  in  the  backscattering  direction  and  gives  rise  to  the  RCS. 
The  backscattered  field  is  calculated  from 


— bs 

E 


e  a0  +  4>’  a* 


dxdy  ejko(ux+vy) 

A 

Ex  cost})1  +  Ey  sine})1 
( -  Ex  sin  cj)1  +  Ey  cos  t})1)  cos  01 


u  =  sin  0‘  cos  (j)1 
v  =  sin  01  sin  (j)1. 
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Ex  and  Ey  are  the  x  and  y  components  of  the  outgoing  field  on  the  aperture  Ea-  The  PO  integral 
defined  above  is  evaluated  by  summing  the  backscattering  contributions  due  to  each  ray  tube. 
The  RCS  is  shown  in  Table  1  [16]. 

Table  1.  Radar  Cross  Section  Definition 


1  Polarization  (1=1, 1=0) 

|  Polarization  (1=0, 1=1) 

Co-pol 

RCS^  =  4k  [A^l2 

RCSee  =  4;t  |A©[2 

Cross-pol 

RCS6(f  =  4k  |Ae|2 

RCS^e  =  4k  \A^\2 

Phase 

Z-A^,  ZA e 

ZAe,  ZA, 

1  =  Amplitude  of  perpendicular  (vertical)  polarization 

I  =  Amplitude  of  parallel  (horizontal)  polarization. 

The  validity  of  XPATCH  results  are  susceptible  to:  (1)  user  selected  input  parameters, 
(2)  measurement  error  of  data  used  in  the  validation  process,  (3)  facetization  levels,  and  (4)  the 
electrical  size  (wavelength  versus  object  size).  All  of  these  factors  are  extremely  important  [9]. 

XPATCH  provides  realistic  RCS,  range  profiles,  and  2-D  SAR  imagery  calculations 
for  complex  targets;  however,  model  deficiencies  exist.  These  include  poor  or  missing  detail, 
inaccurate  surface  curvature,  and  improper  modeling  of  small  features  and  materials.  The  com¬ 
puter  code  does  not  currently  perform  some  of  the  higher  order  scattering  effects  such  as  travel¬ 
ing  waves,  surface  waves,  resonant  effects,  and  creeping  waves  [2].  Upgrading  of  XPATCH  is 
continuous.  Research  is  ongoing  to  quantify  the  level  of  fidelity  required  for  different  frequency 
ranges  and  applications. 

B.  Radar  Cross  Section  -  Basic  Scattering  Code 

ElectroScience  Laboratory  of  the  Ohio  State  University  has  developed  the  user- 
oriented  Radar  Cross  Section  -  Basic  Scattering  Code  (RCS-BSC)  which  provides  far  zone 
scattering  from  complex  structures.  The  first  version  was  published  in  1979,  and  a  second 
version  followed  during  1990,  [24]. 

RCS-BSC  does  not  interface  with  large  CAD  data  descriptions  of  target  geometry. 

All  target  modeling  is  simulated  within  the  program  through  the  use  of  basic  shapes.  The  basic 
shapes  include  multiple-sided  flat  plates,  multiple-section  cone  frustums,  and  finite-composite 
section  ellipsoids  [22].  A  multiple  sided  flat  plate  or  a  box  type  structure  that  is  composed  of  a 
number  of  multiple  sided  plates  attached  together  at  comers  can  be  analyzed  with  all  its  domi¬ 
nant  terms.  All  finite  plates  and  curved  surfaces  have  to  be  perfectly  conducting.  The  infinite 
ground  plane  can  be  dielectric  [17]. 

Once  all  of  the  required  target  geometry  is  input  into  the  code,  RCS-BSC  analyzes  the 
data  and  puts  it  into  the  correct  form  so  that  the  electric  and  magnetic  fields  as  well  as  coupling 
information  can  be  calculated.  A  complete  set  of  incident  and  scattering  polarizations  is  analyzed 
[17].  The  analysis  is  made  from  uniform  asymptotic  techniques  formulated  from  GTD.  The  prin¬ 
cipal  result  is  found  by  summing  the  diffracted  field  from  every  comer  or  edge  visible  to  the  in¬ 
cident  and  observation  directions  [17]  corner  diffraction  coefficients  are  based  on  Michaeli’s 
equivalent  currents  [8].  The  simulation  code  contains  double  diffracted  fields  from  the  corner  to 
edge,  edge  to  comer,  and  comer  to  comer.  It  also  contains  edge  wave  effects  and  multiple  plate 
to  plate  interactions  [17].  Overall,  the  code  encompasses  finding  closed  form  solutions  for  PO 
and  currents  and  asymptotically  matching  them  to  UTD  solutions  and  formats  [17]. 
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The  code  calculates  the  incident  field  as 


Both  incident  polarizations  are  calculated  at  the  same  time.  The  incident  magnitude,  E‘0,  is 
assumed  to  be  one  volt/meter. 

Reference  8  provides  detailed  derivations  of  all  theoretical  techniques  used  in  RCS- 
BSC.  These  include:  (1)  the  corner  diffraction  coefficients;  (2)  the  relation  that  exists  between 
electromagnetic  surface  currents  on  a  wedge  and  equivalent  currents  in  free  space  that  would 
radiate  at  high  frequencies  in  the  same  way  as  the  wedge  surface  currents;  (3)  the  electromagnet¬ 
ic  fields  that  result  when  a  plane  wave  impinges  on  a  wedge;  (4)  Michaeli’s  fringe  and  PO  equiv¬ 
alent  currents  used;  and  (5)  the  comer  diffraction  coefficients  that  result  from  the  use  of  Michae¬ 
li’s  equivalent  currents  to  determine  the  comer  diffraction  coefficients. 

RCS-BSC  was  designed  to  complement  other  techniques.  It  is  a  fast  and  cost  effective 
means  of  optimizing  parameters  and  anticipating  problems  at  the  early  design  stages  of  a  system. 
The  code  can  handle  only  a  limited  number  of  scatterers  and  has  tight  constraints  on  required 
surface  types;  therefore,  it  is  not  considered  the  best  code  for  complex  targets. 


C.  Simulated  Radar  Image 


Another  fully  polarimetric  computer  code  available  today  for  high-frequency  RCS 
prediction  is  the  SRIM  code.  The  code  has  been  under  development  for  several  years  with 
particular  emphasis  being  placed  on  the  simulation  of  polarimetric  sensor  systems.  Reference 
7  has  outlined  the  simulation  process  shown  in  Figure  6. 


Figure  6.  SRIM  Processing 
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The  first  step  is  the  construction  and  sampling  of  a  three-dimensional  model  of  the 
target  and  scene  geometry.  The  target  is  modeled  using  constructive  solid  geometry  and  bound¬ 
ary  representation.  The  ground  region  is  specified  in  terms  of  the  contour  of  its  underlying 
smooth  surface.  Ray  tracing  samples  target  surfaces  and  generates  surface  information  for  the 
scattering  calculation.  Each  surface  of  the  geometry  model  is  assigned  a  set  of  attributes  which 
determines  the  model  used  to  compute  the  scattered  field.  These  attributes  include  whether  the 
surface  is  smooth  or  rough  and  whether  it  is  perfectly  or  finitely  conducting  [7]. 

Scattering  from  smooth  surfaces  is  computed  from  a  generalization  of  the  PO  approxi¬ 
mation.  Scattering  from  rough  surfaces  is  computed  by  drawing  realizations  from  the  first-  and 
second-order  statistical  moments  of  polarization  scattering  matrix.  These  moments  are  obtained 
from  experimental  data  or  analytical  models  [7]. 

RCS  is  calculated  by  coherent  summation  of  the  scattered  fields.  Three  RCS  values  are 
computed:  isolated  target  RCS,  interacting  target  RCS,  and  ground  terrain  RCS  per  unit  area 
[26].  The  backscattered  field  contribution  is  mapped  to  coordinates  in  the  image  plane  and  con¬ 
volved  with  the  system  impulse  response  to  form  the  complex  polarimetric  image  [7]. 

The  following  mathematical  formulation  [7]  describes  the  SAR  simulation  as  the 
predominant  point  of  view.  Letting  a  and  (3  denote  the  receive  and  transmit  polarizations  of 
the  radar  and  u(t)  the  transmitted  waveform,  the  radar  return  signal  is  modeled  as 

vap(0  =  K  f  Yap(r')  u(t  -  2R/c)  ds' 

J  s 

where 

R  =  Radar  Range 
S  =  Illuminated  Surface 
Yap  =  Scalar  Reflectivity  Density. 

For  smooth  surfaces,  this  will  lead  to  an  explicit  formula  for  the  reflectivity  density. 

Defining  a  dyadic  reflectivity  where 

A  j.  A 

Yap  =  ha  •  Y  • 

A 

ha  =  Magnetic  field  polarization  vector  of  receive  antenna 

A 

hp  =  Magnetic  field  polarization  vector  of  transmit  antenna 

you  obtain: 

Y  (f)  =  -  0 . 5  r  x  rixjl  +  RHj  +  fxnx(l+REjxr 

n  =  Surface  normal, 

r  =  Radar  aspect, 

I  =  Identity  dyadic, 

Rh,  Re  =  Fresnel  reflection  dyadics. 
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SRIM  computes  the  contribution  of  the  nth  reflection  point  to  the  p*  image  pixel  as 


g(fp,  rn)  =  K 


a(fp  -  r')  y  (f)  el2nI°'f  ds' 


&sn 


~  a(rp  -  rn)  e1^  Ta 


a(r)  = 
6Sn  = 
= 

fo  = 


Spatial  impulse  response  of  the  radar, 
Surface  patch  illuminated  at  sample  point, 
2jtf0-r', 


(2/K)  i  • 


For  a  smooth  surface,  the  reflectivity  density  y  is  a  continuous  function  over  the  target 
surface.  The  sampled  reflectivity  T  is  determined  by  integrating  over  a  surface  patch  the  product 
of  y  with  a  position-dependent  phase  function.  An  exact  analytical  evaluation  of  the  integral  is 
used  for  planar  surfaces,  and  a  stationary  phase  evaluation  is  used  for  curved  surfaces  [26].  For 
smooth  surfaces: 


rn  =  K  j  Y  (r')  ei2jtT°'(r'_T'")  ds' . 
hsB 

For  rough  surfaces,  it  is  not  possible  to  define  a  continuous  reflectivity  density.  It  is 
assumed  that  optically  rough  surfaces  result  in  noncoherent  diffuse  scattering.  Using  the  central 
limit  theorem,  the  components  of  the  scattered  field  are  modeled  as  complex  circular  Gaussian 
random  variables. 

r  =  fix  f„ 

Fn  =  Magnetic  field  scattering  matrix. 

SRIM  is  a  more  antiquated  code.  It  contains  multiple  interactions  and  calculates 
reflection;  however,  no  diffraction  calculations  are  included. 

D.  MAX-TRACK 

Under  various  Government  contracts,  Georgia  Tech  Research  Institute  (GTRI)  has  de¬ 
veloped  MAX-TRACK,  which  is  essentially  two  programs  interacting  with  one  another.  A  geo¬ 
metric  description  of  the  target  is  entered  using  MAX;  then  the  model  description  from  MAX  is 
used  as  input  into  the  radar  modeling  program  called  TRACK  [6].  TRACK  uses  PO,  GO,  and 
MEC  to  predict  the  returns  from  the  discrete  reflective  targets  [23].  The  output  from  TRACK 
may  be  used  for  post-processing  into  SAR  or  ISAR  imagery,  range  only  output,  azimuth/eleva¬ 
tion  profile,  raster  scan  image,  doppler  spectra,  or  Doppler  Beam  Sharpening  (DBS)  image  [6]. 
A  general  description  of  each  program  follows. 

MAX  was  developed  to  create  and  modify  models  which  TRACK  could  use  to  predict 
targets.  It  stores  all  the  information  needed  to  accurately  describe  target  geometry  and  material 
properties.  These  models  are  a  collection  of  simple  geometric  shapes  which  describe  the  surface 
geometry  of  a  three-dimensional  object.  Elements  currently  supported  include  generalized  ellip¬ 
tical  cone  frusta,  triangular  and  quadrilateral  flat  plates,  generalized  ellipsoids,  dihedrals  and  tri- 
hedrals,  tophats,  hermite  patches,  and  generalized  elliptical  toroids.  Figure  7  illustrates  these 
simple  element  shapes.  These  elements  are  organized  into  a  hierarchical  structure  of  groups  and 
nodes  [6]. 
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Figure  7.  Element  Shapes  Used  by  MAX 

TRACK  predicts  the  RCS.  After  calculating  the  return  from  the  individual  scattering 
centers  on  the  target,  TRACK  can  either  coherently  sum  the  individual  returns  or  retain  the  com¬ 
plex  polarization  scattering  matrix  from  each  scatterer  for  image  post-processing  [6].  The  scat¬ 
tered  returns  from  each  element  of  the  target  model,  except  multiple  scatter  elements,  are  calcu¬ 
lated  with  the  PO  integral.  Multiple  bounce  reflections  which  occur  for  multiple  scatter  elements 
are  treated  with  a  combination  of  GO  and  PO  techniques.  The  scattered  field  from  straight  edges 
is  calculated  with  MEC  [5]. 

These  formulations  are  modified  by  a  diffuse  scattering  algorithm  which  accounts  for 
the  roughness  of  the  individual  scattering  elements.  From  this,  the  diffuse  scatter  from  dielectrics 
and  conductors  is  calculated  [5]. 

MAX-TRACK  is  predominantly  used  for  simple  target  shapes;  however,  modifica¬ 
tions  are  on-going  for  managing  the  more  complex  targets.  Since  the  code  uses  beam  tracing 
rather  than  ray  tracing,  its  computer  run  time  is  much  faster  than  some  of  the  other  codes. 
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V.  FUTURE  DIRECTION 

Deriving  a  radar  signature  model  through  the  use  of  a  three-dimensional  target  drawing 
database  has  long  been  the  dream  of  radar  developers.  This  type  of  model  would  provide  a  cost 
effective  alternative  to  real  data  collections  and  would  allow  us  to  exploit  operational  parame¬ 
ters.  Through  today’s  computational  ability,  the  dream  of  an  automated  and  systematic  approach 
to  theoretical  RCS  prediction  is  possible. 

Future  direction  for  high-frequency  target  modeling  includes  adding  background  in  the 
models,  development  of  a  near-field  version,  and  development  of  a  hybrid  code  where  different 
materials  such  as  the  wings  on  an  airplane  or  blades  on  a  helicopter  can  be  incorporated.  Reduc¬ 
tion  of  computer  run  time  continues  to  be  crucial.  Many  of  the  existing  codes  take  weeks  to  pre¬ 
dict  the  RCS  for  one  complex  target.  Accuracy  is  critical. 

A  great  deal  of  effort  has  been  placed  in  standardizing  and  validating  the  existing  RCS 
simulation  codes.  An  organization  chartered  as  the  Electromagnetic  Code  Consortium  (EMCC) 
was  formed  by  representatives  of  three  U.  S.  DOD  services  and  NASA  in  1987  to  integrate 
U.  S.  efforts  on  the  development  of  electromagnetic  codes,  and  to  be  a  focal  point  for  advancing 
the  state-of-the-art  in  RCS  code  development  and  research  [25].  The  organization’s  efforts  in¬ 
clude  participation  from  industry  as  well  as  university  scientists  and  engineers. 

Modern  high-frequency  codes  are  migrating  toward  the  use  of  ray  tracing  as  a  means  of 
computing  the  scattered  field  from  complex  objects.  Using  this  technique,  XPATCH  is  consid¬ 
ered  by  the  EMCC  as  the  code  to  use  today  for  high-frequency  modeling  [21].  Refinements,  as 
well  as  additional  capabilities,  are  being  made  continuously. 

Extensive  effort  continues  in  the  development  of  new  codes  as  well  as  improving  existing 
codes.  These  models  must  to  be  accurate,  efficient,  and  user  friendly  [20]. 
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